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The conductance of disordered wires with symplectic symmetry is studied by numerical 
simulations on the basis of a tight-binding model on a square lattice consisting of M lattice 
sites in the transverse direction. If the potential range of scatterers is much larger than the 
lattice constant, the number N of conducting channels becomes odd (even) when M is odd 
(even). The average dimensionless conductance (g) is calculated as a function of system length 
L. It is shown that when N is odd, the conductance behaves as (g) — > 1 with increasing L. 
This indicates the absence of Anderson localization. In the even-channel case, the ordinary 
localization behavior arises and (g) decays exponentially with increasing L. It is also shown 
that the decay of (<?) is much faster in the odd-channel case than in the even-channel case. 
These numerical results are in qualitative agreement with existing analytic theories. 
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The concept of universality classes is important in 
considering quantum electron transport in quasi-one- 
dimensional disordered wires. 1 The universality classes 
describe transport properties which are independent of 
the microscopic details of disordered wires. Ordinary dis- 
ordered wires are classified into either of three univer- 
sality classes (orthogonal class, unitary class and sym- 
plectic class), called the standard three, according to 
the symmetries which they possess. The orthogonal class 
consists of systems having both time-reversal symme- 
try and spin-rotation invariance, while the unitary class 
is characterized by the absence of time-reversal symme- 
try. The systems having time-reversal symmetry without 
spin-rotation invariance belong to the symplectic class. 
We here briefly summarize the typical differences among 
the standard three classes. If system length L is much 
shorter than the localization length £, we say that the 
system is in the metallic regime. In this regime, the weak- 
localization (weak-antilocalization) effect arises in the or- 
thogonal (symplectic) class, while the effect disappears 
in the unitary class. If L is much longer than £, the con- 
ductance decays exponentially with increasing L. This 
localization behavior is observed in all the classes. How- 
ever, the localization length £ normalized by the mean 
free path I is shortest (longest) for the orthogonal (sym- 
plectic) class among the standard three. That is, the 
localization behavior is most evident in the orthogonal 
class. 

Exceptional behaviors, which are incompatible with 
the standard three, have been found in the electron trans- 
port properties of metallic carbon nanotubes (CNs). 2 ~ 8 
Ando and co-workers 2,3 studied CNs of length L when 
the potential range of scatterers is larger than the lat- 
tice constant. They proved that one perfectly conduct- 
ing channel is present even when L — ► oo, based on the 
peculiar symmetry of the reflection matrix. This indi- 
cates the absence of Anderson localization even in the 
limit of L — > oo. It should be noted that a graphite sheet 



has symplectic symmetry 9 and that the number of con- 
ducting channels is always odd in CNs with long-range 
scatterers. Thus, the anomalous behavior should be at- 
tributed to the intrinsic nature of the symplectic class 
with an odd number of channels. This indicates that we 
must separately treat the odd-channel case and the or- 
dinary even-channel case in considering the symplectic 
class. To clarify the nature of the symplectic class with 
an odd number of channels, the present author 5 derived 
the so-called DMPK equation for the transmission eigen- 
values based on a random-matrix theory. It is shown that 
the conductance decays exponentially as g — > 1 with in- 
creasing L. It is also shown that the decay is character- 
ized by the typical length scale £ dd = (1/4) (AT — 1)1, 
where N is the number of conducting channels. Noting 
that the localization length for the even-channel case is 
Coven = (2 AT — 2)l, x we observe that the decay of g in the 
odd-channel case is much faster than that in the even- 
channel case. The same result has been obtained by a 
supersymmetric field approach. 8 

The predicted transport properties are so anomalous 
for the symplectic class with an odd-number of channels 
that further studies based on another independent ap- 
proach are needed to reexamine the predictions. In doing 
so, an approach which is applicable to both the odd- and 
even-channel cases is highly desirable because it enables 
us to examine the predicted even-odd differences. In this 
letter, we consider a tight-binding model on a square lat- 
tice, which consists of M lattice sites in the transverse 
direction, as a model for disordered wires with symplectic 
symmetry. We show that if the potential range of scatter- 
ers is much larger than the lattice constant, the number 
A~ of conducting channels becomes odd (even) when M 
is odd (even). Thus, we can consider not only the odd- 
channel case but also the even-channel cases within this 
model. We calculate the ensemble-averaged dimension- 
less conductance (g) as a function of system length L. It 
is shown that when A" is odd, the conductance behaves 
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as (g) — > 1 with increasing L. In the even-channel case, 
the ordinary localization behavior arises and (g) vanishes 
with increasing L. It is also shown that the decay of (g) 
with increasing L is much faster in the odd-channel case 
than in the even-channel case. These numerical results 
are qualitatively in agreement with the DMPK equation. 

As a microscopic model for disordered wires with sym- 
plectic symmetry, we consider the following tight-binding 
model on a square lattice, which is infinitely long in the 
longitudinal direction and consists of M lattice sites in 
the transverse direction, 

oo M 

*= E E E 

j=-oo 1 = 1 cr,CT'=T,i 

x {\j,l + l,a)V x (a,a')(j,l,a'\+h.c.} 

oo M 

+ E E E 

j=—oo 1 = 1 cr,<7'=T,| 

x {\j + 1, 1, a)V v (a, a')(j, I, a'\ + h.c.} 

oo M 

+ E E E \j^^)V^ P (j,l)(j,l,a\ (1) 

j=-oo 1=1 cr=T,i 
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Here, Vi mp (j, I) denotes the impurity potential at the 
(j, l)th site. Note that our Hamiltonian has full symplcc- 
tic symmetry and can be viewed as a discrete version 
of the effective-mass Hamiltonian for a graphite sheet 
if t x = t y . 2 It is equivalent to Ando model 10 in the 
strong spin-orbit interaction limit. We impose the peri- 
odic boundary condition in the transverse direction. That 
is, \j, M + 1, a) = \j, 1, a). Wc hereafter assume that the 
lattice constant is equal to unity. 

As a function of the longitudinal wave number q (— ir < 
q < 7r), the dispersion relation is 



E n (q) = ±2t y Wsin 2 q + r 2 sin 2 



2irn 



(4) 



where r — t x /t y , n = 0, ±1, • • • , ±(M — l)/2 for an odd 
M and n = 0, ±1, • • • , ±(M - 2) /2, M/2 for an even M. 
We hereafter consider only the subbands with E n {q) > 0. 
We see that there are two valleys: one located at q = 
and the other at q = ±7r. Let us call the former one valley 
A and the latter one valley B. We consider the odd-M 
case. The nondegenerate lowest subband with n — pro- 
vides a conducting channel for each valley if the Fermi en- 
ergy Ep satisfies < Ep < 2t y . The other subbands are 
two-fold degenerate. This means that each subband pro- 
vides two conducting channels for each valley if it crosses 
the Fermi level. Thus, the total number 7V t otai of conduct- 
ing channels is equal to 2, 6, 10, 14, 18, • • • depending on 
Ep. In the even-M case, the lowest subband corresponds 



to n = and M/2 and is thus two-fold degenerate. Fur- 
thermore, the other subbands are four-fold degenerate in 
this case. Consequently, iV to tai = 4, 12, 20, • • • . 

The above argument tells us that the total number 
of conducting channels is always even. However, there is 
a simple way to effectively realize the odd-channel case. 
The following way is inspired by the observation that the 
odd-channel case can be realized in CNs if the potential 
range of scatterers is larger than the lattice constant. 2,3 
We focus our attention on the case where Ep is very 
small, that is, < Ep <C 2t y . In this case, q for any con- 
ducting channels in valley A becomes very small, while 
that in valley B is nearly equal to ±7r. This means that in- 
tervalley scattering must be accompanied by a large mo- 
mentum transfer. Such a scattering is suppressed when 
the spatial range of scatterers is much larger than the 
lattice constant. Thus, valley A and valley B effectively 
decouple each other if only such scatterers are present. 
This means that if M is odd, we can realize a disordered 
wire with symplectic symmetry having an odd number of 
conducting channels. Let N be the number of conducting 
channels for each valley. N — 1, 3, 5, 7, 9, • • • depending 
on Ep in the odd-M case, while N = 2, 6, 10, • • • in the 
even-M case. Thus, we can examine the even-odd differ- 
ences within a single model. 

We assume that the impurity potential at the site (j, I) 
arising from a scatterer at the site (jo,^o) is given by 

u -jo? (i-i r 



V( jo ,i )U,l) - w (joM exp 



A 2 

9 



(5) 

To suppress the intervalley scattering, A y must be larger 
than the lattice constant. We randomly distribute scat- 
terers in the finite region of 1 < j < L. Thus, Vi mp (j, I) 
in the Hamiltonian is given by 



L M 

v imp (j,i)= ]T ^2v UoM {j,i). 

jo = l lo=l 



(6) 



The amplitude w/j^ is distributed uniformly within the 
range of — W/2 < wun < W/2. Let p be the probability 
that each site is occupied by a scatterer. The strength 
of the impurity potential is controlled by W and p. We 
consider that the length of our system is equal to L. 

The number of conducting channels in each valley is 
N, so that the total dimensionless conductance is given 

by 



N 



Stotal 



E E i*-, 

a, 5=1 X,Y=A,B 



Y6| , 



(7) 



where tx a ,Yb is the transmission coefficient for an elec- 
tron coming from the 6th channel in valley Y and going 
to the ath channel in valley X. If the two valleys are 
completely decoupled, eq. (7) is reduced to 



JV 



5total = 2 E I^Aa.Afcl 2 - 



(8) 



a,b=l 



Expecting that a good decoupling is achieved, we focus 
on the dimensionless conductance g for each valley, which 
is defined as g = (/total/2. We can numerically obtain 
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txa,Yb by a recursive Green-function technique. We cal- 
culate the ensemble average (g) for the case of M = 17 as 
an example of the odd-channel case and that for M = 16 
as an example of the even-channel case. In the following 
numerical calculations, we set r = 0.1, A y =4, = 1, 
W/ty — 0.7 and p = 0.2. The ensemble average is taken 
over 1000 samples for all the cases. We here present the 
reason why we have set r = 0.1. As noted above, E F 
should satisfy < E F <C 2t y to suppress the interval- 
ley scattering. If M is very large, we can realize a given 
channel number N for a sufficiently small E F . However, 
if M is not large, the condition < E F <C 2t y cannot be 
satisfied for a given N and thereby the intervalley scat- 
tering is relatively enhanced. In such a situation, we can 
suppress it by reducing r. Note that with decreasing r, 
the lowest subband does not change, while the dispersion 
for each higher subband goes down towards that for the 
lowest subband. Thus, a given N can be achieved by a 
relatively small E F compared with the isotropic case of 
r = 1. 



(9) 
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Fig. 1. The dimensionless conductance (g) for N ■ 
a function of system length L. 



5, 7 and 9 as 



The ensemble-averaged dimensionless conductance (g) 
for the odd-channel case is displayed in Fig. 1 as a func- 
tion of L. The Fermi energies are chosen as E F /t y = 
0.08, 0.11 and 0.14, which correspond to N = 5, 7 and 9, 
respectively. We observe that (g) decreases towards unity 
with increasing L. This result is in agreement with the 
DMPK equation which predicts the conductance quan- 
tization (i.e., g = 1) in the long-L limit. The deviation 
(Sg) = (g) — 1 from the quantized value is displayed in 
Fig. 2. We observe that {Sg) decays exponentially with 
increasing L and that the decay becomes faster with de- 
creasing N. For the even-channel case, we calculated the 
averaged dimensionless conductance (g) for E F /t y = 0.08 
and 0.15, which correspond to N = 6 and 10, respec- 
tively, and display the results in Fig. 2. The behavior of 
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Fig. 2. The deviation (Sg) for N = 5, 7 and 9 and the dimension- 
less conductance (g) for N = 6 and 10 as a function of system 
length L. 



the conductance in the even-channel case is very similar 
to that of (Sg). It decays exponentially with increasing 
L and the decay becomes faster with decreasing N. This 
tendency is in agreement with the DMPK equation. 5 

Now we compare the behavior of (Sg) in the odd- 
channel case and that of (g) in the even-channel case. 
From Fig. 2, we observe that the decay of (g) is much 
slower than that of (Sg). This tendency is also in agree- 
ment with the DMPK equation. Let us estimate char- 
acteristic length scales. The L-dependence of (Sg) is as- 
sumed to be 



(Sg) 



(9) 



where we call £ oc id the conductance decay length. For the 
even-channel case, we assume that 



(g) 



~ e 2 «o 



(10) 



where (even is the localization length. The DMPK equa- 
tion predicts 5 



while 1 



UW = \(N-l)l, 



£cvc„(A0 = (2A^ - 2)1, 



(11) 



(12) 



where I is the mean free path. From the above equations, 
we find that 



£cvcn(10) 



= 9. 



(13) 



W9) 

We obtain from our numerical results that £ e ven(10) « 
260 and £ odd (9) w 70, which result in £cvon(10)/£odd(9) « 
3.7. Thus, our numerical results are not in quantitative 
agreement with the DMPK equation, although we have 
obtained the qualitative agreement. The reason for the 
quantitative discrepancy is not clear at this stage, how- 
ever, it may be related to vertex corrections. Suzuura and 
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Ando 9 have calculated the weak-antilocalization correc- 
tion in a graphite sheet by a diagrammatic perturbation 
theory and have shown that vertex corrections are impor- 
tant for quantitative argument. In contrast, the deriva- 
tion of the DMPK equation is based on the assumption 
of equivalent scattering channels, 5, 13 within which such 
corrections vanish. It is very interesting to study the in- 
fluence of vertex corrections on (g) in the long-i regime, 
however, this is difficult. 

In summary, we have considered a tight-binding model 
on a square lattice as a model for disordered wires which 
belong to the symplectic universality class. We have 
shown that if the potential range of scatterers is much 
larger than the lattice constant, this model enables us 
to study both the odd- and even-channel cases. We have 
calculated the ensemble-averaged dimensionless conduc- 
tance (g) as a function of system length L. It is shown 
that the conductance in the odd-channel case behaves 
as (g) — » 1 with increasing L, while (g) — > in the 
even-channel case. It is also shown that the decay of {g) 
with increasing L is much faster in the odd-channel case 



than in the even-channel case. These numerical results 
are qualitatively in agreement with the DMPK equation. 
However, for the typical length scales which character- 
ize the decay of (g) , we have not obtained a quantitative 
agreement between the numerical results and the DMPK 
equation. 
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